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ABSTRACT 


The  time  transformations  and  time  dependent  inputs  for  spatial  trans¬ 
formations  in  the  Advanced  Orbit/Ephemeris  Subsystem  (AOES)  presently  being 
implemented  by  the  Air  Force  Satellite  Control  Facility  have  been  reviewed. 

The  principle  results  are: 

1.  The  relations  described  in  the  Milestone  2  documents  describing 
the  computer  routines  include  all  relations  required  to  satisfy 
the  accuracy  design  goals  of  AOES. 

2.  If  future  systems  require  higher  accuracy  the  present  time-related 
transformations  will  not  be  adequate.  The  improvements  must  involve 
the  following: 

a.  Wander  of  the  pole  of  rotation  with  respect  to  the  crust  of 
the  earth  can  no  longer  be  ignored. 

b.  The  empirical  relations  between  earth  rotation  time  scales 
(UT2,  UT1)  and  atomic  time  scales  (Al,  UTC,  etc.)  must  be 
improved,  probably  by  fitting  over  shorter  periods  of  time. 

c.  Nutation  terms  of  amplitude  less  than  0.2  arc  second  must  be 
included. 

3.  A  best  set  of  constants  to  be  employed  by  the  AOES  relating  UT2  to 
UTC  for  the  period  1  January  1961  to  30  June  1969  have  been 
calculated.  A  procedure  is  suggested  for  the  continuous  review  of 
the  applicability  of  these  and  subsequent  constants,  with  provisions 
for  updating  them  as  required. 

4.  Values  for  comparison  with  computer  routine  results  are  presented 
for  many  of  the  quantities  under  study. 
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I .  INTRODUCTION 


The  spsce-time  reference  frame  most  convenient  for  the  calculations  of 
satellite  ephemerldes  is  not,  because  of  the  earth's  motion,  the  most 
convenient  reference  frame  in  which  satellite  positions  may  be  measured  by 
an  observer  on  the  earth.  Integration  of  the  equations  of  motion  describing 
the  satellite  motion  is  most  conveniently  carried  out  in  an  Inertial 
coordinate  system,  while  observation  is  most  conveniently  carried  out  in  a 
coordinate  system  fixed  with  respect  to  the  surface  of  the  earth. 

Both  inertial  coordinate  systems  and  earth  fixed  systems  are  to  some 
extent  arbitrary  in  their  definitions  of  origins  and  fixed  directions. 

However,  once  these  definitions  are  established,-*  transformation  from  the 
inertial  to  the  earth  fixed  systems  and  vice-versa  is  a  straightforward, 
though  often  tedious,  computational  exercise.  This  transformation  requires 
a  knowledge  of  the  earth's  motion  with  accuracy  commensurate  with  the  accuracy 
required  from  the  transformation. 

The  present  report  will  concentrate^ attention  on  the  coordinate  systems 
and  tranformatlon  schemes  employed  by  the  Advanced  Orblt/Ephemerls  Subsystem 
(AOES)  computer  programs  presently  being  implemented  by  the  Air  Force 
Satellite  Control  Facility.  The  relations  between  the  time  scales  employed 
by  it  and  the  computation  of  time  dependent  inputs  to  the  spatial  transfor¬ 
mation  procedures  of  AOES  will  be  reviewed.  The  basic  equations  will  be 
presented  and  the  numerical  constants  appearing  in  the  programs  will  be 
evaluated.  Where  possible,  test  vslues,  against  which  computer  routines 
may  be  checked,  will  be  presented. 


II.  AGES  SPACE-TIME  COORDINATE  SYSTEMS 


There  are  two  space-time  reference  systems  used  in  AOES.  The  atomic 
tim*.  and  inertial  system  is  employed  for  orbit  calculations  while  the  earth 
fixed  system  with  a  precisely  defined  civil  time  (UTC)  is  employed  in  the 
rest  of  the  Satellite  Control  Facility.  In  this  section  we  discuss  the  basis 
for  each  of  the  two  spatial  reference  frames  used  and  separately  the  basis 
for  the  two  time  systems. 

A.  SPATIAL  COORDINATE  SYSTEMS 

1.  INERTIAL  COORDINATE  SYSTEM 

An  inertial  coordinate  system  may  be  defined  as  one  in  which  Newton's 
second  law  of  motion  (force  equals  rate  of  change  of  momentum)  holds.  The  AOES 
inertial  coordinate  system  is  defined  with  its  origin  at  the  center  of  mass 
of  the  earth  and  the  orientation  of  its  axes  defined  by  the  mean  equator  of 
the  year  1950 .0,  and  the  intersection  of  this  equator  with  the  mean  ecliptic 
of  1950,  i.e.,  the  mean  equinox  of  1950.0. 

Actually  the  "inertial"  reference  frame  with  origin  fixed  in  the  moving 
earth  used  by  AOES  is  only  an  approximation  to  an  inertial  reference  frame. 

It  will  be  a  satisfactory  approximation  so  long  as  the  acceleration  of  the 
satellite  due  to  gravitational  attraction  by  the  earth  is  large  compared 
with  the  acceleration  of  the  origin  of  the  coordinate  system.  This  requirement 
is  satisfied  for  earth  orbits  well  beyond  synchronous  altitudes.  This 
encompasses  the  uses  that  will  be  made  of  AOES.  Therefore  we  follow  the 
convention  of  calling  a  coordinate  system  with  origin  fixed  st  the  center  of 
the  earth  and  with  its  directions  fixed  in  space,  an  inertial  coordinate 
system. 
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2.  EARTH  FIXED  COORDINATE  SYSTEM 

The  AOES  earth  fixed  reference  frame  has  Its  origin  at  the  center  of 
mass  of  the  earth  and  its  orientation  specified  by  the  equator  and  prime 
meridian  of  the  earth.  The  AOES  actually  uses  several  different  coordinate 
systems  which  are  related  to  this  one  by  non-time-dependent  relations;  for 
example:  position  determined  as  azimuth,  elevation,  and  range  from  a  known 
spot  on  the  surface  of  the  earth.  Because  we  are  at  present  emphasizing 
only  time  dependent  relations,  all  of  these  earth  fixed  coordinate  systems 
will  be  regarded  as  equivalent. 

B.  TIME  SCALES 

"Time  is  not  physically  tangible;  it  has  no  unique  physical  property 
that  permits  its  laboratory  examination.  Time  is  esaentially  metaphysical; 
there  is  no  direct  way  of  measuring  it,  even  in  principle.  Nonetheless,  our 
lives  are  ordered  by  it  and,  more  Importantly,  our  physics  is  also  ordered 
by  it.  Therefore  time  must  be  measured.  The  foundations  of  ph,  leal  science 
include  the  article  of  faith  that  a  'uniform*  time  exists  that  corresponds 
Identically  with  the  variable  called  time  in  dynamics.  Even  Isaac  Newton 
suspected  (Ref.  1,  page  8)  the  impossibility  of  determining  'uniform*  time, 
and  he  commented  upon  the  necessity  of  distinguishing  between  this  construct 
of  faith  and  the  physical  measure?  of  time.  The  failure  to  reconcile 
observations  with  dynamical  theories  may  lead  one  to  amend  or  discard  the 
laws  of  dynamics  or  the  means  of  determination  of  time,  but  the  faith  in  a 
'uniform*  supertime  is  untouched"  (Ref.  2,  page  14).  In  the  following  we 
review  the  specific  definitions  of  the  two  time  scales  employed  in  AOES. 

1.  ATOMIC  TIME,  A1 

Because  the  calculation  of  satellite  ephemerldcs  is  based  on  this  faith 
in  the  existence  of  a  "uniform"  supertime,  it  is  convenient  to  use  a  practical 
time  scale  which  agrees  as  far  as  is  known  with  the  uniform  time.  The  A1 
atomic  time  scale  is  the  one  presently  available  that  meets  this  requirement 
and  has  been  chosen  as  the  time  scale  for  the  inertial  space-time  coordinate 
system  in  AOES.  Ue  briefly  review  the  development  of  this  time  scale. 
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The  positions  of  astronomical  objects  have  been,  until  recently,  the 
j  Is  for  the  practical  determination  of  uniform  time.  These  positions 
have  been  employed  because  they  could  be  measured  quite  precisely  and 
because  classical  dynamics,  which  was  the  first  branch  of  physics  to  develop 
e  precise  mathematical  formalism,  was  adequate  to  predict  most  of  these 
positions  with  high  precision.  Since  the  chief  cause  for  apparent  motion 
of  astronomical  objects  Is  tb.?  motion  of  the  earth,  these  time  measuring 
schemes  where  In  fact  based  on  the  motion  of  the  earth. 

As  more  precise  independent  measurements  of  time  became  possible  by 
means  of  quarts  oscillators  and  atomic  beam  systems,  Irregularities  were 
discovered  In  time  based  on  earth  motion  which  could  not  be  quantitatively 
predicted  but  could  be  qualitatively  described  and  were  attributed  to  such 
things  as  tidal  motions,  and  seasonal  changes  In  the  polar  Ice  masses.  It 
thus  became  necessary  to  define  tbe  unit  of  time,  the  ephemerls  second,  as 
e  certain  fraction  of  a  particular  year  near  the  beginning  of  this  century 
(Refs.  3,  4).  The  epbemerls  time  defined  in  this  way  can  in  prludple  always 
be  found  from  the  reduction  of  a  sufficient  number  of  star  position  observa¬ 
tions  but  the  determination  is  tedious  and  can  be  done  only  after  an  event  has 
happened. 

The  development  of  the  cesium  beam  clock  In  the  early  1950's  introduced 
a  means  of  keeping  time  more  precise  than  ephemerls  time  and  with  greater 
convenience.  The  ephemerls  second  was  measured  to  be  9  192  631  770  periods 

of  the  microwave  trensltlon  between  tbe  hyperflne  levels  F  »  4,  ■  0  and 

2  133 

F  •  3,  Dp  ■  C  of  the  Sj^  ground  state  of  the  unperturbed  Cs  atom 

(Ref.  5).  A  new  time  scale  (Al)  has  been  defined  with  the  above  number  of 

oscillations  defined  as  one  second  and  0b  0m  0s  of  1  January  1958  defined 
h  n  s 

as  0  0  0  of  the  etomlc  time  scale.  This  causes  Al  to  differ  from  ephemerls 

* 

time  by  32.15  second.  Up  to  tbe  present  tbls  difference  has  remained 


A 

The  best  and  most  recent  velue  of  ET  -  Al  Is  32.15  seconds  (Ref.  17). 

Earlier  meesurements,  however,  gave  a  value  of  32.25  seconds.  This  value 
la  still  used  in  the  calculation  of  many  Jet  Propulsion  Laboratory  ephemerldes 
(Ref.  2).  Since  the  only  use  of  ephemerls  time  In  AOES  Is  believed  to  be 
obtaining  quantities  from  these  JPL  ephemerldes  It  Is  desirable  tbat  the 
same  constent  be  used  to  relate  £  and  Al.  This  requires  that  ETTOAT,  the 
only  an try  In  CBLK,  'AURLK2,  be  32.25  seconds.  If  JPL  changes  the  value 
used  in  their  ephemerls  computations  then  ETTOAT  should  be  changed. 
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constant  indicating  that,  within  the  precision  of  measurements  made  to  date, 

ephemeris  time  and  atomic  time  are  equivalent.  Recognizing  this  and  the 

convenience  with  which  time  can  be  determined  from  atomic  standards  the 

second  was  redefined  in  1967  as  9  192  631  770  periods  of  the  above 
133 

mentioned  Cs  resonance  line.  It  is  this  atomic  time  (Al)  which  is  used 
in  the  AOES  as  the  best  approximation  to  uniform  time  and  thus  the  time 
variable  of  integration  for  determining  satellite  position. 

2.  UNIVERSAL  TIME 

Although  atomic  time  is  the  most  convenient  for  orbital  calculations,  a 
time  scale  related  to  the  rotation  of  the  earth  is  employed  by  the  Air  Force 
for  observations  and  communication  with  other  Satellite  Control  Facility 
Systems. 

The  time,  using  the  earth  as  the  clock,  may  be  determined  from  the 
observed  position  of  some  known  point  among  the  stars.  More  specifically  the 
angle  along  the  celestial  equator  from  the  great  circle  through  the  observer's 
zenith  west  to  the  great  circle  through  the  known  point  is  used.  These 
angles,  usually  measured  in  units  of  hours,  are  called  hour  angles.  The 
local  time,  XT,  on  the  x  scale,  where  x  is  the  accepted  reference  point  in 
the  sky,  is  then  the  hour  angle  of  reference  point;  XT  «  ^a(x)*  By  th*s 
definition  only  observers  at  the  same  longitude  will  have  the  same  time, 
which  is  inconvenient.  Therefore  a  single  meridian  is  defined  by  convention 
as  the  time  keeping  meridian.  Then  the  local  time,  XT^^  ■  XT  +  longitude, 
where  XT  is  the  time  at  the  prime  meridian  and  longitude  is  measured  increasing 
to  the  west  from  the  prime  meridian. 

Two  reference  points  and  their  related  time  scales  are  required  in  AOES. 
One  point  is  the  equinox  of  date  and  the  associated  time  scale  is  called 
sidereal  time  (ST) =>  This  time  (or  angle)  is  required  to  relate  the  "inertial" 
coordinate  system  to  the  earth  fixed  one.  The  second  reference  point,  which 
defines  the  universal  time  scale,  the  basis  for  civil  time  keeping,  is  always 
very  nearly  directly  opposite  the  sun.  Because  of  the  earth's  motion  about 


the  sun,  the  second  reference  point  is  in  constant  motion  with  respect  to  the 
first.  A  relation  between  them  is  required.  By  definition: 


CT  ~  U- 

““(equiapx) 

The  angle  from  the  equinox  to  the  wiversal  time  point  is  defined  (Ref.  6) 

to  be  12  hours  +  R  where 
u 

R  -  18h  38™  45*836  +  8640184! 542  T  +  o!o929T2 
u 


and  T  is  the  number  of  Julian  centuries  of  36525  days  which  have  elapsed 
since  12h  UT  on  0  January  1900. 


UT  «  ha  «  ha  -  R  -  12n 
u  eq  u 

h  (1) 

UT  -  ST  -  R  -  12 
u 

UT  is  thus  determined  by  means  of  Eq.  1  from  sidereal  time  obtained  by  observing 
the  zenith  transit  of  stars  whose  position  with  respect  to  the  equinox  is 
accurately  known. 

When  this  time  scale  is  compared  with  an  independent  precise  time  scale, 
such  as  atomic  time,  UT  is  observed  to  have  a  varying  difference  from  the 
atomic  time.  There  appears  to  be  1)  changes  which  depend  on  the  observer's 
location  and  result  from  the  wandering  of  the  earth's  pole  of  rotation  with 
respect  to  the  crust  of  the  earth;  2)  changes  which  are  seasonal;  and 
3)  gradual  slowing  down  of  the  earth  rotation  rate.  None  of  these  can  be 
completely  predicted  on  a  theoretical  basis  and  are  the  subject  of  study  due 
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to  their  intrinsic  interest  (Refs.  7,  8).  Three  time  scales  In  use  reflect  these 
various  differences:  UTU  is  the  time  scale  which  results  frou  applying  Eq.  1 
to  neasured  star  positions.  UT1  is  the  result  of  correcting  UTO  for  polar 
notion,  UT2  is  the  result  of  correcting  UT1  for  seasonal  variations.  The  only 
effect  left  in  UT2,  then  is  the  gradual  slowing  down  of  the  earth. 

The  UTC  tine  scale  is  designed  to  provide  a  convenient  tine  scale  with 
the  stability  of  atonic  clocks  and  still  be  a  close  approximation  to  UT2. 

Because  of  the  measurements  and  corrections  required,  UT2  can  be  determined 
only  long  after  an  event,  and  is  therefore  inconvenient  for  many  purposes. 

The  ratio  of  the  length  of  the  UTC  second  to  the  A1  second  length  is  held 
constant  at  a  value  such  that  the  UTC  second  is  approximately  the  same  length 
as  the  UT2  second.  The  value  of  this  ratio  is  decided  upon  and  announced 
in  advance  by  an  international  bureau  on  the  basis  of  recent  trends  in  UT2. 

When  the  difference  between  UT2  and  UTC  becomes  larger  than  about 

100  milliseconds,  a  step  is  put  in  UTC  to  bring  th'.  two  closer  together. 

UTC  is  the  time  scale  used  by  AOES  for  communication  with  other  Satellite 
Control  Facility  systems. 

It  is  important  to  recognize  that  UTC  and  Ai  have  a  mathematically  defined 
relation  between  them.  Given  one  the  other  can  always  be  found  to  within 
the  accuracy  of  constants  in  a  formula.  Similarly  UTO  and  sidereal  time  have 
a  mathematically  defined  relationship.  On  the  other  hand  the  relations 
between  UTO  UT1,  UT2,  and  UTC  are  basically  empirical,  though  they  can  be 
estimated  with  fair  accuracy  on  the  basis  of  past  experience. 


III.  COORDINATE  SYSTEM  TRANSFORMATIONS 


Relating  the  position  of  a  object  as  expressed  in  inertial  coordinates 
to  the  position  in  earth  fixed  coordinates  requires  only  a  knowledge  of  the 
orientation  of  one  coordinate  system  with  respect  to  the  other  at  any  given 
time.  Since  the  largest  part  of  the  change  in  orientation  is  due  to  the 
earth's  motion,  which  is  for  the  most  part  quite  predictable,  the  most 
logical  way  to  obtain  the  required  transformation  would  seem  to  be  to 
compute  the  relative  orientation  at  any  specific  time,  to  v^rly  this  to  relate 
the  position  vector  in  the  two  coordinate  systems,  and  finally  to  apply  small 
corrections,  predicted  as  best  as  possible  by  statistical  means  on  the  basis 
of  recent  empirical  data  on  the  motion  of  the  earth. 

This  is  generally  the  procedure  actually  followed  in  AOES  although  the 
detailed  way  in  which  the  corrections  are  made  reflects  the  traditional 
coordinate  systems  of  the  astronomer.  The  transformation  is  broken  into 
several  parts,  some  of  which  (precession  and  nutation)  are  introduced  as 
changes  in  the  inertial  coordinate  system,  others  (rotation  of  earth  about 
its  axis,  seasonal  and  random  variations  in  the  rotation  rate)  are  introduced 
via  the  time,  and  some,  such  as  polar  wander  are  applied  directly.  In  addition 
to  the  spatial  transformations  there  is  a  time  transformation  between  A1  and 
UTC. 

The  steps  in  the  transformations  between  the  inertial  system  of  AOES  and 
the  earth  fixed  systems  of  AOES  are  indicated  schematically  in  Fig.  1.  The 
mathematical  formulation  of  the  transformations  is  described  in  detail  in 
Ref.  9  and  will  not  be  repeated  in  its  entirety  here.  The  purpose  of  the 
present  report  is  1)  to  review  the  inputs  to  these  transformations  which  are 
time  dependent,  and  2)  to  review  the  transformations  between  the  A1  time 
scale  used  in  the  inertial  system  and  the  UTC  time  scale  used  in  communicating 
with  the  world  outside  the  AOES.  These  are  the  operations  taking  place  within 
the  broken  line  in  Fig.  1. 
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Schematic  Representation  of  Space~Time  Reference  Systems 
in  AQES  Emphasizing  Transformations  Between  Them  which 
are  Time  Dependent.  Those  Computations  Referred  to 
within  tl  e  Broken  Line  are  the  Detailed  Concern  of  this 
Report. 


The  AOES  computer  routines  which  perform  these  calculations  are:  ’DATU 
(Atomic  to  Universal  Time),  ’DOME  (Sidereal  Time),  and  'DPNUT  (Precession 
Nutation  Matrix).  A  fourth  routine,  'DATXME,  performs  housekeeping  functions 
such  as  converting  the  time  between  two  dates  into  minutes.  Since  It 
Introduces  no  new  time  scales  or  transformations.  It  will  not  be  reviewed 
here.  The  physical  constants  required  for  ’DATU  and  'DOME  are  found  In  the 
system  data  base  as  a  constant  block  ’AUBLK.  In  the  following  sections 
numerical  values  of  all  the  constants  In  'AUBIK  will  be  calculated.  This 
Information  is  summarized  In  Appendix  D,  which,  for  every  entry  in  ’AUBLK, 
either  lists  the  value  calculated  In  this  report  or  refers  to  the  tabulation 
of  that  value  elsewhere  In  the  report.  The  constants  required  for  ’DPNUT  are 
part  of  the  code  for  that  routine. 

The  accuracy  of  the  tranformatlon  should  meet  and  exceed  the  design  goals 
of  AOES.  In  terms  of  a  nominal  90  minute  orbit  "nothing  in  the  computer 
program  shall  preclude  the  attainment  of  these  accuracies:"  50  ft  cross 
track,  200  ft  In  track  and  100  ft  radial  distance  (Ref.  10).  In  angular 
measure  the  50  ft  requirement  for  a  satellite  at  185  miles  with  a  90  minute 
period  requires  an  angular  accuracy  of  0.47  arc  sec  or  30  milliseconds  time. 

A.  TIME  DEPENDENT  SPATIAL  TRANSFORMATION  INPUTS 

Referring  to  Fig.  1,  there  are  two  time  dependent  Inputs  to  the  spatial 
transformations.  These  are  the  sidereal  time  (computed  by  ’DOME)  and  the 
precession-nutation  matrix  (computed  by  ’DPNUT) . 

1.  SIDEREAL  ANGLE 

The  most  important  quantity  required  for  transformation  between  the  earth 
fixed  and  the  AOES  Inertial  system  Is  the  angle  between  the  prime  meridian 
of  the  earth  fixed  system  and  the  meridian  through  the  mean  equinox  of  date. 
This,  however,  is  just  the  hour  angle  of  the  equinox,  or  the  sidereal  time. 
From  Eq.  1  of  the  earlier  discussion  of  time  bases,  this  is: 

e  -  ST  »  UT  +■  R  +  12h 
u 
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Converting  to  degrees  and  AOES  base  time,  tbis  becomes: 

0  -  UT  +  100!0755416  +  0t9856473458  d  +  2.90  x  10-13d2  (2) 

* 

wbere  d  is  the  universal  time  in  days  since  AOES  reference  time.  Some  checks 
of  the  computer  routine  for  sidereal  angle  have  been  conducted  (Ref.  11). 
Values  for  additional  checks  may  be  obtained  from  tbe  ephemeris  (Ref.  12)  for 
any  year  under  the  heading  Mean  Sidereal  Time. 

2.  PRECESSION  AND  NUTATION 

Because  of  the  gravitational  attraction  of  other  bodies  in  the  solar 
system  tbe  rotation  axis  of  the  earth  does  not  remain  fixed  in  space,  but 
slowly,  over  a  period  of  about  26000  years,  describes  s  circle  of  about 
23.5*  radius.  Furthermore  the  combined  effects  of  the  planets  serve  to  slowly 
change  the  plane  of  tbe  earth's  orbit  (tbe  ecliptic).  The  combined  motion 
is  known  as  general  precession  of  the  mean  pole.  In  sddition,  there  are 
vsrious  shorter  period  motions  of  the  true  pole  about  the  m^n  pole  position 
described  as  nutation. 

Ue  are  principally  concerned  here  that  the  angular  arguments  usea  to 
compute  the  nutation  and  precession  matrices  are  correct.  In  order  to  be 
specific  about  the  srguments,  however,  it  is  necessary  to  give  the  way  in 
which  they  are  used,  which  means  specifying  the  transformation  matrices. 

Assume  we  wish  to  transform  a  vector  specified  by  its  position  with  respect 
to  the  mean  equinox  and  true  equator  of  date,  (center  of  Fig.  1)  to  a  posi¬ 
tion  vector  in  the  AOES  inertial  reference  frame  which  referred  to  the 
mean  equinox  and  equator  of  1950  (left  side  of  Fig.  1).  Transformation  from 


To  avoid  very  large  numbers  in  some  calculations  all  times  in  AOES  are,  by 
convention,  measured  from  a  reference  time  of  0^  1  January  1950. 

The  first  two  constants  In  Eq.  2,  expressed  as  fractions  of  a  revolution, 
l’JKl  aid  THK2  in  constant  block  ’AUBLK  in  the  AOES  data  bsse.  For  some 
purposes  the  rate  of  change  of  sidereal  time  with  respect  to  UT2  is  needed. 
This  may  be  obtained  by  differentiating  the  above  equation.  The  result  is 
THDKX  in  'AUBLK.  In  all  cases  the  quadratic  term  is  ignored. 
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mean  equinox  and  true  equator  of  date  to  mean  equinox  and  mean  equator  of 
date  is  accomplished  by  the  nutation  matrix  N.  Transformation  from  mean 
equinox  and  equator  of  date  to  mean  equinox  and  equator  of  1950.0  is  then 
accomplished  by  the  precession  matrix,  P. 


-*•  -*• 
r  -  PN  r 
o 

We  look  at  each  matrix  in  some  detail  now. 

a.  Nutation 

Nutation  is  specified  by  two  angles:  Aiji,  the  angular  separation  measured 
along  the  ecliptic  between  the  mean  and  true  equinoxes  (K  and  K*  in  Fig.  2), 
and  Ae,  the  difference  in  angle  between  the  equatorial  planes  and  *-tae 
ecliptic  plane  for  the  true  equator  e^  and  the  me^n  equator  The  trans¬ 
formation  is  first  Atjf  cos  e  about  z,  then  e^  about  the  new  x  axis,  then  Aiji 
about  the  latest  z  axis,  and  finally  (-e^)  about  the  latest  x  axis.  Combining- 
all  these  operations  into  one  matrix,  using  the  small  angle  approxi&iation  and 
keeping  terms  only  to  the  first  order  in  Ac  and  Ai|>,  the  resulting  nutation 
matrix  N  is 

/  X 

N  - 1  0  1 

\-A4»  sin  E  -Ae 

The  arguments  required  are  AtJ>  and  Ae  which  are  given  by  trigonometric  series, 
the  constants  of  which  have  been  determined  partially  theoretically  and 
partially  empirically. 

N 

Ai|<  »  l  a.  (t)  sin  [a.  (t)  ] 
i-1  1 
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z 


Figure  2.  Geometry  Relating  the  True  Equator  and  Mean 

Equator  of  Date  from  which  the  Nutation  Matrix 
may  be  Computed.  The  Angle  KK'  is  the  Nutation 
in  Longitude,  A>{>.  The  Difference  Between  the 
Inclination  of  the  True  Equator  to  the  Ecliptic 
(e.)  and  the  Mean  Equator  to  the  Ecliptic  (e.) 
la  the  Nutation  in  Obliquity,  Ac. 
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(13) 


N 

Ae  -  l  b  (t)  cos  (a.Ct)] 

1-1  1  1 

A 

where  t  is  time.  The  a  ■  £  (L.  F  (t),  are  various  linear  combinations  of 

1  j-1  J 

Fj(t)  which  are  in  turn  polynomials  in  t.  Listed  in  the  Explanatory  Notes  to 

Ephemerls  (Ref.  6)  are  69  terms  with  a^  0.0002  arc  sec  and  AO  with  b^ 

0.0002  arc  sec.  The  0.A7  arc  sec  accuracy  required  by  A0ES  suggests  that  not 

nearly  all  of  these  terms  are  required.  Table  la  lists  the  coefficients  of 

the  functions  which  have  either  ai  or  bi  greater  than  0.1  arc  sec.  Table  lb 

lists  the  coefficients  in  the  F(t)  polynomials.  The  time  arguments  (d  and 

n 

Tq)  to  be  used  with  the  coefficients  of  Table  1  are  measured  from  the  beginning 

of  the  present  Julian  Century  while  the  time  arguments  in  A0ES  are  measured 

in  days,  d,  since  A0ES  base  time.  The  relations  between  these  quantities 

are:  d  -  d  +  18262.5  and  T  -  d  /36525. 
n  n  n 

Computer  calculations  may  be  checked  by  comparison  with  "he  value.'  listed 

in  Table  2.  For  10  day  intervals  during  1968,  the  nutation  in  longitude  (A>J>) 

and  nutation  in  obliquity  (Ac)  are  listed  as  obtained  both  from  the  limited 

series  of  Table  1  and  from  the  full  series  as  tablulated  in  the  1968  ephemerls 

(Ref.  12).  The  nontrivial  matrix  elements,  N^  ■  Ai{)  sin  e,  and  “  Ac 

corresponding  to  the  limited  series  calculation  are  also  listed.  The  maximum 

difference  between  the  ephemerls  and  the  calculated  values  is  0.171  arc  sec. 

The  average  difference  without  regard  to  sign,  is  0.059  arc  sec,  well  within 

_  * 

the  accuracy  requirement  of  A0ES. 


* 

DPNUT  does  not  include  the  i  ■  A  term  from  Table  la.  When  this  term  is  not 
included,  the  same  comparison  over  1968  gives  a  maximum  difference  of  0.261 
arc  sec  and  an  average  without  regard  to  sign  of  0.083  arc  sec.  The  maximum 
error  is  within  a  factor  of  two  of  the  A0ES  design  limit. 


The  Data  Dynamics  Inc.  literature  on  'DPNUT  chose  to  test  the  nutation  matrix 
on  1.^0  January  1980.  The  values  calculated  on  the  basis  of  Table  1  for  this 
date  are  Aif;  -  -7’.'958,  N._  -  -1.5350  x  ur5,  and  N  _  -  -A.25A7  x  10-5.  If  i  -  A 
term  of  Table  la  is  notincluded,  we  then  calculate  Aip  -  7'.‘952,  Ae  -  -8V776. 

N. _  *  -1.5338  x  10“5,  and  N23  ■  -A.25A7  x  lo-5.  These  latter  values  are  in 
good  agreement  with  the  test  number  published  in  the  8/15/07  Milestone  A, 
revision  B,  of  'DPNUT. 
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b.  Coefficients  of  nutation  polynomials 


Fj  “  fj0  +  fjldn  +  fJ2dn  +  fj3dn 


Index,  J 

fJ0 

(deg) 

fJl 

(deg) 

fJ2 

(deg) 

fJ3 

(deg) 

1 

358.475833 

0.9856002669 

-1.12  x  10-13 

-6.8  x  10-20 

2 

11.250889 

13.2293504490 

-2.407  x  10-12 

-7.0  x  10-21 

3 

350.737486 

12.1907491914 

-1.076  x  10*12 

3.9  x  10‘20 

4 

259.183275 

-0.0529539222 

1.557  x  10-12 

4.6  x  10_2° 
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Table  2 .  Nutation  Matrix  Values  for  1968 


Julian  Day  Nutation  in  Umgitiult  Nutation  in  Obliquity  Off-Axis  Matrix  Elements 

(arc  mc)  (arc  acc)  (Units  of  10  ) 


Calculated 

Epbaesrls 
(Kaf.  13) 

Calculated 

Epheeerls 
(Kef.  13) 

"1.3 

“2.3 

243S855.5 

-6.477 

-6.401 

7.736 

7.708 

-1.2494 

3.7507 

65.6 

-6.149 

-6.132 

7.946 

7.949 

-1.1861 

3.8521 

75.5 

-5.457 

-5.428 

8.203 

8.224 

-1.0527 

3.9768 

85.5 

-4.949 

-4.822 

8.299 

8.320 

-0.9546 

4.0236 

95.5 

-5.027 

-5.089 

8.457 

8.472 

-0.9697 

4.1001 

2439905.5 

-5.098 

-5.150 

8.791 

8.832 

-0.9835 

4.2622 

15.5 

-4.901 

-4.891 

9.008 

9.046 

-0.9454 

4.3674 

25.5 

-5.058 

-5.209 

9.014 

8.996 

-0.9757 

4.3699 

35.5 

-5.630 

-5.683 

9.085 

9.064 

-1.0859 

4.4044 

45.5 

-5.890 

-5.948 

9.227 

9.223 

-1.1361 

4.4734 

2439955.5 

-5.839 

-5.927 

9.150 

9.107 

-1.1264 

4.4359 

65.5 

-6.096 

-5.992 

8.918 

8.874 

-1.1760 

4.3233 

75.5 

-6.423 

-6.444 

8.826 

8.797 

-1.2390 

4.2789 

85.5 

-6.178 

-6.182 

8.776 

8.775 

-1.1918 

4.2549 

95.5 

-5.696 

-5.525 

8.538 

8.543 

-1.0988 

4.1392 

2440005.5 

-5.532 

-5.542 

8.307 

8.292 

-1.0671 

4.0272 

15.5 

-5.249 

-5.292 

8.316 

8.345 

-1.0125 

4.0316 

25.5 

-4.450 

-4.362 

8.355 

8.389 

-0.8584 

4.0506 

35.5 

-3.745 

-3.793 

8.267 

8.261 

-0.7225 

4.0079 

45.5 

-3.496 

-3.542 

8.311 

8.303 

-0.6744 

4.0294 

2440055.5 

-3.117 

-3.119 

8.585 

8.592 

-0.6013 

4.1619 

65.5 

-2.478 

-2.525 

8.788 

8.782 

-0.4780 

4.2604 

75.5 

-2.271 

-2.213 

8.861 

8.833 

-0.4381 

4.2958 

85.5 

-2.522 

-2.533 

9.075 

9.046 

-0.4865 

4.3996 

95.5 

-2.54o 

-2.567 

9.383 

9.383 

-0.4912 

4.5490 

2440105.5 

-2.452 

-2.329 

9.478 

9.477 

-0.4730 

4.5953 

15.5 

-2.842 

-2.835 

9.437 

9.407 

-0.5483 

4.5754 

25.5 

-3.391 

-3.378 

9.520 

9.516 

-0.6540 

4.6153 

35.5 

-3.464 

-3.375 

9.575 

9.601 

-0.6683 

4.6419 

45.5 

-3.426 

-3.436 

9.376 

9.375 

-0.6608 

4.5458 

2440155.5 

-3.723 

-3.649 

9.142 

9.142 

-0.7183 

4.4323 

65.6 

-3.811 

-3.804 

9.078 

9.103 

-0.7352 

4.4013 

75.5 

-3.309 

-3.370 

8.961 

8.989 

-0.6383 

4.3446 

85.5 

-2.827 

-2.729 

8.681 

8.701 

-0.5454 

4.2086 

95.5 

-2.627 

-2.667 

8.531 

8.529 

-0.5067 

4.1362 

2440205.5 

-2.088 

-2.196 

8.598 

8.629 

-0.4028 

4.1686 

15.5 

-1.173 

-1.092 

8.599 

8.624 

-0.2263 

4.1689 

25.5 

-0.609 

8.530 

-0.1175 

4.1354 
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b.  Precession 

The  precession  is  described  by  the  three  angles  shown  in  Fig.  3.  a  »90*-£o 
is  the  right  ascension  of  the  ascending  node  (n)  of  the  equator  of  epoch  t 
on  the  equator  of  tg  reckoned  from  the  equinox  of  tQ.  8  ■  90°  +  z  is  the 
right  ascension  of  the  node  (n)  reckoned  from  the  equinox  of  t  on  the 
equator  of  t  and  6  is  the  inclination  of  the  equator  of  t  to  the  equator  of  tg. 
To  transform  X,  Y,  Z  components  of  one  system  to  XQ,  Yg,  and  Zg  of  the  other 
requires  first  a  rotation  +  8  about  Z,  -0  about  the  new  X,  and  finally  -a  about 
the  latest  Z  axis.  The  P  matrix  elements  then  result  from  straightforward 
trigonometry. 

P..  ■  cos  S  cos  6  cos  Z  -  sin  £  sin  Z 
11  o  o 

P21  =  -sin  £o  cos  6  cos  Z  -  cos  5  sin  Z 

P^  ■  -sin  0  cos  Z 

P,_  *  cos  5  cos  0  sin  Z  +  sin  i;  cos  Z 
1/  o  o 

P00  ■  -sin  Zy  cos  0  sin  Z  +  cos  5  cos  Z 
Li  O  O 

Pj2  *  -sin  0  sin  Z 

P^  ■  cos  ?o  sin  6 

P2j  ■  -  sin  sin  0 

P33  *  cos  ® 

The  angles  are  given  by  (Ref.  6): 

C  -  (2304V250  +  1V396  T_  )  T  +  0V  302T2  +  0'.'018  T3 
0  Op  p  p  p 

z  -  C  +  O'.1 791  T2 
o  p 

0  «  (2004 V682  -  01853  T_  )  T  -  O'. '426  T2  -  0  V 042  T3 

Op  p  p  p 
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Tp  is  the  time  interval  between  the  two  epochs  measured  in  tropical  centuries 
of  36524.22  days.  TQp  is  the  time  since  Besselian  year  1900.0,  also  measured 
in  tropical  centuries.  Changing  to  AOES  base  time,  which  is  0.077  day^fnto 
Bessellian  year  1950,  TQp  -  0.5  and  Tp  -  (0.077  +  d)/36524.22,  where  d  is  the 
number  of  days  since  AOES  reference  time.  The  equations  then  become: 

£  “  2304V948  T  +  0V302  T2  +  0V018  T3 
o  P  P  P 

z  -  £  +  O'. '791  T2 
o  p 

6  -  2004'. '2555  T  -  0V426  T2  -  0V042  T3 
P  P  P 

The  results  of  calculating  the  angles  and  the  matrix  elements  for  a  few 
selected  times  are  tabulated  in  Table  3.  Also  tabulated  are  results  for 
approximately  the  same  times  taken  from  Tables  2.1  and  2.2  of  the 
Explanatory  Notes  to  the  Ephemeris  (Ref.  6).  The  differences  in  the  angles 
are  due  to  the  calculation  being  made  for  0^0  on  January  1  while  the  ephemeris 

•Jf 

tables  were  computed  for  the  beginning  of  the  Besselian  Year. 

B.  TIKE  TRANSFORMATIONS 

The  time  transformations  required  are  between  UTC,  used  in  earth  fixed 
system,  and  A1  used  in  the  inertial  system.  Since  sidereal  time  is  directly 
proportional  to  CTO,  the  relation  between  UT0  and  A1  or  UTC  must  also  be 
accurately  known.  Universal  time  is  required  to  calculate  the  precession- 
nutation  matrix  but  since  this  is  a  small  correction  the  differences  between 
UTC  and  A1  are  insijtinif icant .  This  section  reviews  these  transformations 
between  time  scales. 


The  times  for  which  values  are  preseuted  in  Table  3  are  the  times  used  for 
testing  the  P  matrix  in  'DPNUT.  The  values  in  Table  3  for  the  off  axis 
matrix  elements  differ  slightly  but  significantly  from  the  values  published 
in  'DPNUT,  Milestone  4,  Revision  B,  8/15/67.  It  is  suspected  this  is  due  to 
roundoff  error  in  the  Milestone  4  values  but  this  is  not  definitely  known. 


Table  3.  Calculated  Freceaelon  Matrix  Valuee 


The  UTC  time  scale  is  related  to  the  A1  scale  as  follows: 

A1  -  UTC  +  a£  -  Si  (UTC  -  UTC^ 

The  s^  are  frequency  offsets  decided  in  advance  and  announced  by  an  inter¬ 
national  bureau  to  keep  the  rate  of  change  of  UTC  with  respect  to  A1  about 
equal  to  the  average  rate  of  change  of  UT2  with  respect  to  Al.  The  are 
changed  only  on  the  beginning  of  a  year.  The  a^  are  also  decided  upon  and 
announced  in  advance  by  the  same  international  bureau,  when  UTC  and  UT2  differ 
by  more  than  about  0.1  sec.  The  a^  may  be  changed  at  the  beginning  of  any 
month. 

Values  of  the  frequency  shifts,  Sj^  and  offsets,  a^,  have  been  calculated 
for  all  the  changes  which  have  taken  place  in  UTC  between  1  January  1961  and 
30  June  1969.  These  are  presented  in  the  first  three  columns  of  Table  4. 

The  values  of  UTC^^  appearing  in  the  first  column  are  given  in  units  of  days 
since  AOES  base  time.  (The  1  August  1963  entry  is  an  exception  in  that  no 
defined  change  in  UTC  took  place  at  this  time.  This  entry  was  introduced  to 
reduce  the  maximum  difference  between  UTC  and  UT2  (See  Appendix  C) . 

During  the  early  part  of  this  time  period  the  clocks  used  could  not  be 
held  at  precisely  the  nominal  announced  frequency  offset,  This  is 

reflected  in  Table  4  by  the  s^  which  do  not  have  a  simple  value.  The  values 
for  s^  which  are  specified  to  3  decimal  place?  are  the  results  of  least  squares 
fitting  to  the  time  scale  established  by  the  National  Bureau  of  Standard  Radio 
Station  (WW)  broadcasts.  The  UTC  values  calculated  from  the  constants  differ 
in  all  cases  by  Less  than  one  millisecond  from  the  UTC  time  scale  established 
by  WWV.  The  advances  in  frequency  control  at  U.S.  Naval  Observatory  (USNO) 
and  within  the  Satellite  Control  Facility  should  make  corrections  to  the 
nominal  offset,  s^,  unnecessary  in  the  future. 


Table  4.  AUTAB  from  Least  Squares  Fit  to  US  Naval  Observatory  Data 
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Timing  Polynomials,  A1  to  UT1 

As  discussed  under  rime  scales,  the  relationship  between  atomic  basjd 

time  scales,  such  as  UTC  or  Ai,  and  earth  motion  time  scales  such  as  UTO  is, 

when  considered  at  their  highest  accuracy,  a  relation  which  is  known  precisely 

only  after  an  event  has  occurred.  A  most  straightforward  way  to  obtain  UTO 

from  Al  would  be  a  power  series  fit  to  obtain  UT1,  then  the  station  peculiar 

correction  to  get  UTO.  The  constants  in  the  formula  would  be  changed 

whenever  the  fit  became  unsatisfactory.  This  is  the  approach  taken  by  the 

Jet  Propulsion  Laboratory  (Refs.  14,  IS)  and  the  coefficients  of  their 

timing  polynomials  are  readily  available  (Table  3,  Ref.  2).  From  these 

coefficients  "the  value  of  Al-UTl  is  known  to  ±0.005  sec.  after  thft  fact 

2 

and  to  0.0002  sec/day  additional  when  predicting  ahead."  To  accomplish  this 
accuracy  requires  41  different  sets  of  polynomials  to  cover  the  period  from 
1  January  1961  through  June  1969. 

b.  Successive  Corrections,  UTC  to  UT2 

Another  approach  is  to  take  advantage  of  AOES's  use  of  UTC,  which  is 
already  a  first  approximation  to  UT2  and  apply  in  reverse  the  corrections 
used  to  get  UTC  from  UT1.  This  is  the  approach  followed  by  AOES.  A  linear 
relation  between  UTC  and  UT2  is  assumed. 

UT2  -  UTC  +  C±  +  d±  (UTC  -  UTC^ 

As  a  programming  convenience  the  same  UTC^  used  for  the  UTC  to  Al  conversion 
are  used  here.  Tabulated  in  Table  4  are  the  constants  resulting  from  a  least 
squares  fit  of  this  formula  to  the  USNO  determined  UT2  scale.  (Ref.  16).  The 
maximum  values  of  |UT2(calculated)  -  UT2(USNC)|  are  also  tabulated  in  Table  4- 
the  absolute  difference  between  the  UT2  calculated  from  the  constants  and  the 
UT2  established  by  the  US  Naval  Observatory. 


The  largest  errors  are  around  20  msec,  some  four  times  larger  than  the 
JPL  timing  polynomials  errors,  but  only  14  different  sets  of  constants  are 
required  to  cover  the  1961  to  1969  period  compared  vith  the  41  required  for 
JFL  polynomial .  The  22  msec  is  still  well  within  the  30  msec  accuracy 
requirement  of  AOES.  In  Appendix  A,  all  of  the  USNO  time  data  employed 
are  tabulated  along  with  the  calculated  values  of  UTC  and  UT2  obtained  from 
the  constants  of  Table  4  for  each  input  time. 

c.  UT2  to  UT1 

After  determining  UT2,  the  seasonal  correction,  S,  must  be  applied  to 
obtain  UT1,  UT1  •  UT2  -  S.  The  assumed  form  of  S  is: 

S  ■  Uj^  sin  (2nt)  +  u^  cos  (2irt)  +  u^  sin  (4irt)  +  u^  cos  (4itt) 

where  t  is  the  fraction  of  a  Julian  year  of  365.25  days  since  1  January  of 

* 

the  year.  The  u^  are  chosen  for  an  entire  year.  Values  found  by  a  least 

squares  fit  to  the  USNO  data  (Ref.  16)  are  listed  in  Table  5.  The  complete 

USNO  data  and  the  fit  for  all  input  values  is  shown  in  Appendix  B.  For  all 

years  except  1962  the  errors  are  less  tban  the  1  msec  accuracy  with  which 

the  USNO  values  are  stated. 


Table  5.  Seasonal  Time  Variation  Coefficients 


Year 

U1 

(sec) 

1961 

0.0221 

1962 

0.0210 

1963-68 

0.0220 

*The  u  are  tbe  U1K1,  U1K2,  UJK3 
in  the  AOES  data  base. 


U2 

U4 

(sec) 

(aec) 

(sec) 

-0.0169 

-0.0069 

0.0059 

-0.0135 

-0.G074 

0.0067 

-0.0120 

-0.0059 

0.0070 

U1K4  of  the 

constant  block 

'AUBLK 
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d.  UT1  to  iiTO 

The  correction  from  UT1  to  UTO  .'t>  really  another  coordinate  change  which 
should  be  applied  to  relate  earth  fixed  coordinates  referred  to  a  point  on  the 
earth  to  earth  fixed  coordinates  referred  to  the  pole  of  rotation.  If  r  is  a 
vector  determined  with  respect  to  a  system  of  coordinates  defined  by  the  mean 
pole  of  1903  and  r'  is  a  vector  determined  with  respect  to  a  system  of 
coordinates  defined  by  the  current  pole  of  rotation  then 

-+■ 

r  ■  R  r 

where 


And  x  and  y  are  shifts  from  the  mean  pole  in  radian  measure.  If  A  is  the 
longitude  and  0  the  latitude,  this  same  change  can  be  represented  as  a  charge 
of  latitude,  A<t>  and  longitude,  AA. 

AA  *  A '-A  -  tan  4>(x  sin  A  +  y  cos  A) 

A<p  -  <$'-<()  “  x  cos  A  -  y  sin  A 

The  values  of  x  and  y  are  presented  as  a  polar  plot  for  the  1962  to  1968 
time  period  in  Fig.  4  (Ref.  2).  It  will  be  noted  that  over  a  period  of 
6  months  the  values  have  changed  by  as  much  as  0.5  arc  sec.  This  is 
approximately  equal  to  the  design  goal  of  AOES  and  is  not  corrected  for  in 
AOES  tranformations  since  the  correction  is  station  peculiar. 
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e  Accuracy 

The  largeat  single  error  source  is  the  polar  motion  which  is  not 

corrected  for  and  may  amount  to  ±.020  sec  time.  The  next  largest  source 

of  error  ia  the  converalon  from  UTC  to  UT2  which  may  at  times  be  in  error 

nearly  as  much.  The  other  corrections  are  never  in  error  by  more  than  0.1 

of  these  amounts.  If  these  two  sources  of  error  were  Independent  and  random 

(neither  of  which  is  probably  true)  then  one  might  expect  the  total  error 

1/2 

to  be  (2)  times  the  error  from  either  of  the  two  sources  of  equal  magni¬ 
tude  or  about  ±  0.030  sec  time — just  about  the  A0ES  design  goal. 

Thus  the  present  AOES  transformation  achieves  the  design  accuracy 
goals.  Improvement  in  accuracy  of  other  parts  of  the  system  beyond  the 
present  goals  is  not  useful  unless  the  transformations  are  also  improved. 

For  greater  accuracy  the  relation  between  A1  (or  UTC)  and  UT2  must  be  improved, 
perhaps  by  uaing  something  like  the  JPL  timing  polynomials,  and  the  polar 
motion  must  be  included. 

3.  PARAMETER  UPDATE 

Moat  of  the  parameters  occurring  in  the  various  time  transformation 
equations  are  not  varying,  and  except  to  review  them  occasionally  to  assure 
that  the  best  currently  available  values  are  in  use,  they  do  not  need  to 
be  changed.  Exceptions  to  this  are  the  constants  relating  A1  to  UTC  to  UT1. 
Theae  constants  must  b_  updated  to  reflect  the  most  recent  available 
information.  A  convenient  source  of  this  current  information  is  the  weekly 
bulletin  of  the  US  Naval  Observatory,  "Preliminary  Times  and  Coordinates  of 
the  Pole,  Series  7",  which  gives  information  on  UTC,  UT2,  UT1,  UT0  and  should 
it  be  useful  sometime  in  the  future  x,  and  y  coordinates  of  the  pole. 

Appendix  C  outlines  a  possible  systematic  way  to  incorporate  this  data  as  it 
becomes  available,  while  at  the  same  time  making  a  minimum  number  of  changes 
in  the  computer  data  base. 
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x  (arc  sec) 

POSITIVE  TO  0°,  NEGATIVE  TO  I80°E ,  180°  W 


y  (arc  sec) 

POSITIVE  TO  90° W, NEGATIVE  TO  90° E 


Figure  A.  Polar  Orbit  Plot  of  the  Coordinates  of  the  Pole  of  Rotation 
of  the  Earth  for  the  Period  1962  to  1968  (Ref.  2). 
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APPENDIX  A 


Comparison  of  UTC  and  UT2  as  computed  from  the  constants  of  Table  4  as 
compared  with  the  US  Naval  Observatory  established  scales,  which  are  labelled 
"EXP". 
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APT ESDI!  3 

Comparison  of  seasonal  correction  S  *  PT2  -  CT1  as  obtained  froa  the 
parameters  of  Table  5  with  the  correction  established  by  the  US  Naval 
Observatory  and  indicated  by  “EXP". 
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APPENDIX  C 


Time  conversion  constant  update  procedure. 

All  of  the  parameters  required  to  convert  from  A1  to  UTC,  from  UTC  to  UT2 
and  from  UT2  to  UT1  are  included  In  a  constant  block  (CBLK)  called  'AUBLK  in 
the  AOES  System  II  Data  Base.  The  constants  for  the  first  two  conversions  are 
contained  In  Table  'AUTAB  within  this  block.  The  correlation  between  the 
symbols  used  in  this  report  and  the  item  names  in  this  table  are  indicated 
by  the  headings  in  Table  4.  'AUTAB  has  space  for  a  maximum  of  10  entries. 

Table  4  already  has  14  entries  so  not  all  of  Table  4  can  be  used  as  'AUTAB. 
Because  of  the  way  'AUTAB  is  employed  by  the  computer  routine  'DATU,  the  entries 
must  occur  in  consecutive  order  with  no  omissions.  Thus  to  support  currently 
flying  missions  the  last  ten  entries  in  Table  4  would  be  used. 

New  entries  must  be  made  in  'AUTAB  whenever  a  step  change  is  made  in  UTC 
or  whenever  the  frequency  offset  Is  changed.  It  may  be  advantageous  to  make 
additional  entries  when  it  is  impossible  to  fit  UTC  to  UT2  within  the  required 
accuracy  over  the  time  period  between  UTC  changes  with  a  single  set  of 
Constanta.  In  the  eight  year  period  covered  by  Table  4  this  has  been 
necessary  only  once,  namely  the  entry  of  1  August  1963. 

In  addition  to  adding  new  entries  to  'AUTAB  the  values  in  the  last  (most 
recent)  entry  should  reflect  as  large  an  amount  of  data  as  is  available  on 
the  relation  between  UTC,  UT2.  The  most  convenient  source  of  this  information 
is  the  "Preliminary  Times  and  Coordinates  of  the  Pole.  Series  7"  published  by 
the  US  Naval  Observatory.  The  method  by  which  this  information  is  Incorporated 
in  the  base  involves  tradeoffs  between  frequent  changes  in  the  data  base, 
level  of  management  authority  which  should  permit  changes  in  the  data  base, 
auxiliary  computer  programs  required,  etc.,  all  considerations  outside  the 
scope  of  this  report.  We  suggest  one  scheme  simply  to  indicate  the  type  of 
procedure  which  might  be  used. 


This  updating  procedure  is  outline  aa  a  flow  chart  in  Fig.  C-l.  The 
criteria  for  making  a  change  in  'AUTAB  is  that  for  some  input  day  there  is 
a  greater  than  30  msec  difference  between  UT2  as  calculated  from  the  constants 
in  'AUTAB  and  UT2  as  obtained  from  the  Naval  Observatory.  If  a  change  ia 
necessary  the  first  attempt  is  to  change  the  constanta  relating  UTC  to  UT2 
to  better  values.  If  elimination  of  the  troublesome  time  difference  ia  not 
possible  this  way  then  a  new  entry  should  be  made  in  'AUTAB. 

A  new  entry  is  made  in  'AUTAB  whenever  the  ]UT2cal~UT2yg{}0  j  ia  too 
large  or  whenever  a  UTC  discontinuity  occurs.  For  either  cause  the  new 
entry  is  always  the  last  one  in  the  table  and  if  one  must  be  dropped  to  keep 
within  the  table  size  it  must  be  the  first  entry  in  the  table.  The  items 
in  the  new  entry  are  determined  as  shown  in  Table  C-l.  When  a  new  entry  is 
made  there  will  be  no  further  need  to  update  the  old  last  entry.  Therefore 
the  best  values  of  UCTU2C  and  UCTU2D  for  that  time  period  should  be  inserted 
even  if  the  old  values  were  giving  acceptable  answers. 


C-2 


MAKE  NEW  ENTRY 
IN  'AUTAB  WITH 
FIRST  DAY  OF 
PRECEEOING  MONTH 
AS  UTCOIS  (COLUMN 
H,  TABLE  C-2) 


MAKE  NEW  ENTRY 
IN  'AUTAB  WITH 
EFFECTIVE  OATE 
OF  UTC  SHIFT  AS 
UTCDIS  (COLUMN 
I  .TABLE  C -I ) 


NO  FURTHER 
ACTION  REQUIRED 


Figure  C-l.  Flow  Chart  of  a  Possible  Procedure  for 
Updating  ’AUTAB. 
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APPENDIX  D 


Susnary  of  Constant  Block  AUBLK  Data  Valuea. 


Name 

Description 

Value  or  Source 

AUTA.I 

UTC  to  A1  and  A1  to  UTC  Conversion 

Table  4  and  Appendix  C 

UTCIND 

Internal  Flag 

A1IND 

Internal  Flag 

U1K1 

U1K2 

\  Seasonal  Variation  Constants 

Table  5 

U1K3  | 

U1K4 

THK1 

Constant  for  Sidereal  Time  Comp. 

0.2779876155  j 

THK2 

Constant  for  Sidereal  Time  Comp. 

-2  ( 

0.2737909294  x  10  c  ;2q. 

THDK1 

Sidereal  Rate  Conatant 

1.002737909294  ) 
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